3.2.100 \(\int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx\) [200]

3.2.100.1 Optimal result
3.2.100.2 Mathematica [A] (verified)
3.2.100.3 Rubi [A] (verified)
3.2.100.4 Maple [C] (verified)
3.2.100.5 Fricas [A] (verification not implemented)
3.2.100.6 Sympy [F]
3.2.100.7 Maxima [A] (verification not implemented)
3.2.100.8 Giac [F(-2)]
3.2.100.9 Mupad [F(-1)]

3.2.100.1 Optimal result

Integrand size = 22, antiderivative size = 160 \[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=\frac {x \sqrt {c+a^2 c x^2}}{24 a^3}-\frac {x^3 \sqrt {c+a^2 c x^2}}{20 a}-\frac {2 \sqrt {c+a^2 c x^2} \arctan (a x)}{15 a^4}+\frac {x^2 \sqrt {c+a^2 c x^2} \arctan (a x)}{15 a^2}+\frac {1}{5} x^4 \sqrt {c+a^2 c x^2} \arctan (a x)+\frac {11 \sqrt {c} \text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {c+a^2 c x^2}}\right )}{120 a^4} \]

output
11/120*arctanh(a*x*c^(1/2)/(a^2*c*x^2+c)^(1/2))*c^(1/2)/a^4+1/24*x*(a^2*c* 
x^2+c)^(1/2)/a^3-1/20*x^3*(a^2*c*x^2+c)^(1/2)/a-2/15*arctan(a*x)*(a^2*c*x^ 
2+c)^(1/2)/a^4+1/15*x^2*arctan(a*x)*(a^2*c*x^2+c)^(1/2)/a^2+1/5*x^4*arctan 
(a*x)*(a^2*c*x^2+c)^(1/2)
 
3.2.100.2 Mathematica [A] (verified)

Time = 0.14 (sec) , antiderivative size = 105, normalized size of antiderivative = 0.66 \[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=\frac {a x \left (5-6 a^2 x^2\right ) \sqrt {c+a^2 c x^2}+8 \sqrt {c+a^2 c x^2} \left (-2+a^2 x^2+3 a^4 x^4\right ) \arctan (a x)+11 \sqrt {c} \log \left (a c x+\sqrt {c} \sqrt {c+a^2 c x^2}\right )}{120 a^4} \]

input
Integrate[x^3*Sqrt[c + a^2*c*x^2]*ArcTan[a*x],x]
 
output
(a*x*(5 - 6*a^2*x^2)*Sqrt[c + a^2*c*x^2] + 8*Sqrt[c + a^2*c*x^2]*(-2 + a^2 
*x^2 + 3*a^4*x^4)*ArcTan[a*x] + 11*Sqrt[c]*Log[a*c*x + Sqrt[c]*Sqrt[c + a^ 
2*c*x^2]])/(120*a^4)
 
3.2.100.3 Rubi [A] (verified)

Time = 0.76 (sec) , antiderivative size = 299, normalized size of antiderivative = 1.87, number of steps used = 13, number of rules used = 12, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.545, Rules used = {5481, 262, 262, 224, 219, 5487, 262, 224, 219, 5465, 224, 219}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x^3 \arctan (a x) \sqrt {a^2 c x^2+c} \, dx\)

\(\Big \downarrow \) 5481

\(\displaystyle \frac {1}{5} c \int \frac {x^3 \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx-\frac {1}{5} a c \int \frac {x^4}{\sqrt {a^2 c x^2+c}}dx+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}\)

\(\Big \downarrow \) 262

\(\displaystyle \frac {1}{5} c \int \frac {x^3 \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \int \frac {x^2}{\sqrt {a^2 c x^2+c}}dx}{4 a^2}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}\)

\(\Big \downarrow \) 262

\(\displaystyle \frac {1}{5} c \int \frac {x^3 \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\int \frac {1}{\sqrt {a^2 c x^2+c}}dx}{2 a^2}\right )}{4 a^2}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}\)

\(\Big \downarrow \) 224

\(\displaystyle \frac {1}{5} c \int \frac {x^3 \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\int \frac {1}{1-\frac {a^2 c x^2}{a^2 c x^2+c}}d\frac {x}{\sqrt {a^2 c x^2+c}}}{2 a^2}\right )}{4 a^2}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {1}{5} c \int \frac {x^3 \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 5487

\(\displaystyle \frac {1}{5} c \left (-\frac {2 \int \frac {x \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx}{3 a^2}-\frac {\int \frac {x^2}{\sqrt {a^2 c x^2+c}}dx}{3 a}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 262

\(\displaystyle \frac {1}{5} c \left (-\frac {2 \int \frac {x \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx}{3 a^2}-\frac {\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\int \frac {1}{\sqrt {a^2 c x^2+c}}dx}{2 a^2}}{3 a}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 224

\(\displaystyle \frac {1}{5} c \left (-\frac {2 \int \frac {x \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx}{3 a^2}-\frac {\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\int \frac {1}{1-\frac {a^2 c x^2}{a^2 c x^2+c}}d\frac {x}{\sqrt {a^2 c x^2+c}}}{2 a^2}}{3 a}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {1}{5} c \left (-\frac {2 \int \frac {x \arctan (a x)}{\sqrt {a^2 c x^2+c}}dx}{3 a^2}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}-\frac {\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}}{3 a}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 5465

\(\displaystyle \frac {1}{5} c \left (-\frac {2 \left (\frac {\arctan (a x) \sqrt {a^2 c x^2+c}}{a^2 c}-\frac {\int \frac {1}{\sqrt {a^2 c x^2+c}}dx}{a}\right )}{3 a^2}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}-\frac {\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}}{3 a}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 224

\(\displaystyle \frac {1}{5} c \left (-\frac {2 \left (\frac {\arctan (a x) \sqrt {a^2 c x^2+c}}{a^2 c}-\frac {\int \frac {1}{1-\frac {a^2 c x^2}{a^2 c x^2+c}}d\frac {x}{\sqrt {a^2 c x^2+c}}}{a}\right )}{3 a^2}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}-\frac {\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}}{3 a}\right )+\frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {1}{5} x^4 \arctan (a x) \sqrt {a^2 c x^2+c}+\frac {1}{5} c \left (-\frac {2 \left (\frac {\arctan (a x) \sqrt {a^2 c x^2+c}}{a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{a^2 \sqrt {c}}\right )}{3 a^2}+\frac {x^2 \arctan (a x) \sqrt {a^2 c x^2+c}}{3 a^2 c}-\frac {\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}}{3 a}\right )-\frac {1}{5} a c \left (\frac {x^3 \sqrt {a^2 c x^2+c}}{4 a^2 c}-\frac {3 \left (\frac {x \sqrt {a^2 c x^2+c}}{2 a^2 c}-\frac {\text {arctanh}\left (\frac {a \sqrt {c} x}{\sqrt {a^2 c x^2+c}}\right )}{2 a^3 \sqrt {c}}\right )}{4 a^2}\right )\)

input
Int[x^3*Sqrt[c + a^2*c*x^2]*ArcTan[a*x],x]
 
output
(x^4*Sqrt[c + a^2*c*x^2]*ArcTan[a*x])/5 - (a*c*((x^3*Sqrt[c + a^2*c*x^2])/ 
(4*a^2*c) - (3*((x*Sqrt[c + a^2*c*x^2])/(2*a^2*c) - ArcTanh[(a*Sqrt[c]*x)/ 
Sqrt[c + a^2*c*x^2]]/(2*a^3*Sqrt[c])))/(4*a^2)))/5 + (c*((x^2*Sqrt[c + a^2 
*c*x^2]*ArcTan[a*x])/(3*a^2*c) - ((x*Sqrt[c + a^2*c*x^2])/(2*a^2*c) - ArcT 
anh[(a*Sqrt[c]*x)/Sqrt[c + a^2*c*x^2]]/(2*a^3*Sqrt[c]))/(3*a) - (2*((Sqrt[ 
c + a^2*c*x^2]*ArcTan[a*x])/(a^2*c) - ArcTanh[(a*Sqrt[c]*x)/Sqrt[c + a^2*c 
*x^2]]/(a^2*Sqrt[c])))/(3*a^2)))/5
 

3.2.100.3.1 Defintions of rubi rules used

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 224
Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Subst[Int[1/(1 - b*x^2), x], 
x, x/Sqrt[a + b*x^2]] /; FreeQ[{a, b}, x] &&  !GtQ[a, 0]
 

rule 262
Int[((c_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^2)^(p_), x_Symbol] :> Simp[c*(c*x) 
^(m - 1)*((a + b*x^2)^(p + 1)/(b*(m + 2*p + 1))), x] - Simp[a*c^2*((m - 1)/ 
(b*(m + 2*p + 1)))   Int[(c*x)^(m - 2)*(a + b*x^2)^p, x], x] /; FreeQ[{a, b 
, c, p}, x] && GtQ[m, 2 - 1] && NeQ[m + 2*p + 1, 0] && IntBinomialQ[a, b, c 
, 2, m, p, x]
 

rule 5465
Int[((a_.) + ArcTan[(c_.)*(x_)]*(b_.))^(p_.)*(x_)*((d_) + (e_.)*(x_)^2)^(q_ 
.), x_Symbol] :> Simp[(d + e*x^2)^(q + 1)*((a + b*ArcTan[c*x])^p/(2*e*(q + 
1))), x] - Simp[b*(p/(2*c*(q + 1)))   Int[(d + e*x^2)^q*(a + b*ArcTan[c*x]) 
^(p - 1), x], x] /; FreeQ[{a, b, c, d, e, q}, x] && EqQ[e, c^2*d] && GtQ[p, 
 0] && NeQ[q, -1]
 

rule 5481
Int[((a_.) + ArcTan[(c_.)*(x_)]*(b_.))*((f_.)*(x_))^(m_)*Sqrt[(d_) + (e_.)* 
(x_)^2], x_Symbol] :> Simp[(f*x)^(m + 1)*Sqrt[d + e*x^2]*((a + b*ArcTan[c*x 
])/(f*(m + 2))), x] + (Simp[d/(m + 2)   Int[(f*x)^m*((a + b*ArcTan[c*x])/Sq 
rt[d + e*x^2]), x], x] - Simp[b*c*(d/(f*(m + 2)))   Int[(f*x)^(m + 1)/Sqrt[ 
d + e*x^2], x], x]) /; FreeQ[{a, b, c, d, e, f, m}, x] && EqQ[e, c^2*d] && 
NeQ[m, -2]
 

rule 5487
Int[(((a_.) + ArcTan[(c_.)*(x_)]*(b_.))^(p_.)*((f_.)*(x_))^(m_))/Sqrt[(d_) 
+ (e_.)*(x_)^2], x_Symbol] :> Simp[f*(f*x)^(m - 1)*Sqrt[d + e*x^2]*((a + b* 
ArcTan[c*x])^p/(c^2*d*m)), x] + (-Simp[b*f*(p/(c*m))   Int[(f*x)^(m - 1)*(( 
a + b*ArcTan[c*x])^(p - 1)/Sqrt[d + e*x^2]), x], x] - Simp[f^2*((m - 1)/(c^ 
2*m))   Int[(f*x)^(m - 2)*((a + b*ArcTan[c*x])^p/Sqrt[d + e*x^2]), x], x]) 
/; FreeQ[{a, b, c, d, e, f}, x] && EqQ[e, c^2*d] && GtQ[p, 0] && GtQ[m, 1]
 
3.2.100.4 Maple [C] (verified)

Result contains complex when optimal does not.

Time = 0.52 (sec) , antiderivative size = 176, normalized size of antiderivative = 1.10

method result size
default \(\frac {\sqrt {c \left (a x -i\right ) \left (a x +i\right )}\, \left (24 \arctan \left (a x \right ) a^{4} x^{4}-6 a^{3} x^{3}+8 a^{2} \arctan \left (a x \right ) x^{2}+5 a x -16 \arctan \left (a x \right )\right )}{120 a^{4}}-\frac {11 \sqrt {c \left (a x -i\right ) \left (a x +i\right )}\, \ln \left (\frac {i a x +1}{\sqrt {a^{2} x^{2}+1}}-i\right )}{120 a^{4} \sqrt {a^{2} x^{2}+1}}+\frac {11 \sqrt {c \left (a x -i\right ) \left (a x +i\right )}\, \ln \left (\frac {i a x +1}{\sqrt {a^{2} x^{2}+1}}+i\right )}{120 a^{4} \sqrt {a^{2} x^{2}+1}}\) \(176\)

input
int(x^3*arctan(a*x)*(a^2*c*x^2+c)^(1/2),x,method=_RETURNVERBOSE)
 
output
1/120/a^4*(c*(a*x-I)*(I+a*x))^(1/2)*(24*arctan(a*x)*a^4*x^4-6*a^3*x^3+8*a^ 
2*arctan(a*x)*x^2+5*a*x-16*arctan(a*x))-11/120/a^4*(c*(a*x-I)*(I+a*x))^(1/ 
2)*ln((1+I*a*x)/(a^2*x^2+1)^(1/2)-I)/(a^2*x^2+1)^(1/2)+11/120/a^4*(c*(a*x- 
I)*(I+a*x))^(1/2)*ln((1+I*a*x)/(a^2*x^2+1)^(1/2)+I)/(a^2*x^2+1)^(1/2)
 
3.2.100.5 Fricas [A] (verification not implemented)

Time = 0.27 (sec) , antiderivative size = 94, normalized size of antiderivative = 0.59 \[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=-\frac {2 \, {\left (6 \, a^{3} x^{3} - 5 \, a x - 8 \, {\left (3 \, a^{4} x^{4} + a^{2} x^{2} - 2\right )} \arctan \left (a x\right )\right )} \sqrt {a^{2} c x^{2} + c} - 11 \, \sqrt {c} \log \left (-2 \, a^{2} c x^{2} - 2 \, \sqrt {a^{2} c x^{2} + c} a \sqrt {c} x - c\right )}{240 \, a^{4}} \]

input
integrate(x^3*arctan(a*x)*(a^2*c*x^2+c)^(1/2),x, algorithm="fricas")
 
output
-1/240*(2*(6*a^3*x^3 - 5*a*x - 8*(3*a^4*x^4 + a^2*x^2 - 2)*arctan(a*x))*sq 
rt(a^2*c*x^2 + c) - 11*sqrt(c)*log(-2*a^2*c*x^2 - 2*sqrt(a^2*c*x^2 + c)*a* 
sqrt(c)*x - c))/a^4
 
3.2.100.6 Sympy [F]

\[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=\int x^{3} \sqrt {c \left (a^{2} x^{2} + 1\right )} \operatorname {atan}{\left (a x \right )}\, dx \]

input
integrate(x**3*atan(a*x)*(a**2*c*x**2+c)**(1/2),x)
 
output
Integral(x**3*sqrt(c*(a**2*x**2 + 1))*atan(a*x), x)
 
3.2.100.7 Maxima [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.79 \[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=-\frac {1}{120} \, {\left (a {\left (\frac {3 \, {\left (\frac {2 \, {\left (a^{2} x^{2} + 1\right )}^{\frac {3}{2}} x}{a^{2}} - \frac {\sqrt {a^{2} x^{2} + 1} x}{a^{2}} - \frac {\operatorname {arsinh}\left (a x\right )}{a^{3}}\right )}}{a^{2}} - \frac {8 \, {\left (\sqrt {a^{2} x^{2} + 1} x + \frac {\operatorname {arsinh}\left (a x\right )}{a}\right )}}{a^{4}}\right )} - 8 \, {\left (\frac {3 \, {\left (a^{2} x^{2} + 1\right )}^{\frac {3}{2}} x^{2}}{a^{2}} - \frac {2 \, {\left (a^{2} x^{2} + 1\right )}^{\frac {3}{2}}}{a^{4}}\right )} \arctan \left (a x\right )\right )} \sqrt {c} \]

input
integrate(x^3*arctan(a*x)*(a^2*c*x^2+c)^(1/2),x, algorithm="maxima")
 
output
-1/120*(a*(3*(2*(a^2*x^2 + 1)^(3/2)*x/a^2 - sqrt(a^2*x^2 + 1)*x/a^2 - arcs 
inh(a*x)/a^3)/a^2 - 8*(sqrt(a^2*x^2 + 1)*x + arcsinh(a*x)/a)/a^4) - 8*(3*( 
a^2*x^2 + 1)^(3/2)*x^2/a^2 - 2*(a^2*x^2 + 1)^(3/2)/a^4)*arctan(a*x))*sqrt( 
c)
 
3.2.100.8 Giac [F(-2)]

Exception generated. \[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=\text {Exception raised: TypeError} \]

input
integrate(x^3*arctan(a*x)*(a^2*c*x^2+c)^(1/2),x, algorithm="giac")
 
output
Exception raised: TypeError >> an error occurred running a Giac command:IN 
PUT:sage2:=int(sage0,sageVARx):;OUTPUT:sym2poly/r2sym(const gen & e,const 
index_m & i,const vecteur & l) Error: Bad Argument Value
 
3.2.100.9 Mupad [F(-1)]

Timed out. \[ \int x^3 \sqrt {c+a^2 c x^2} \arctan (a x) \, dx=\int x^3\,\mathrm {atan}\left (a\,x\right )\,\sqrt {c\,a^2\,x^2+c} \,d x \]

input
int(x^3*atan(a*x)*(c + a^2*c*x^2)^(1/2),x)
 
output
int(x^3*atan(a*x)*(c + a^2*c*x^2)^(1/2), x)